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Abstract 

In linear anisotropic elasticity, the elastic properties of a medium are described by 
the fourth rank elasticity tensor C. The decomposition of C into a partially symmetric 
tensor M and a partially antisymmetric tensors N is often used in the literature. An 
alternative, less well-known decomposition, into the completely symmetric part S of C 
plus the reminder A, turns out to be irreducible under the 3-dimensional general linear 
group. We show that the SA-decomposition is unique, irreducible, and preserves the 
symmetries of the elasticity tensor. The M/V-decomposition fails to have these desirable 
properties and is such inferior from a physical point of view. Various applications of 
the SA-decomposition are discussed: the Cauchy relations (vanishing of A), the non- 
existence of elastic null Lagrangians, the decomposition of the elastic energy and of the 
acoustic wave propagation. The acoustic or Christoffel tensor is split in a Cauchy and 
a non-Cauchy part. The Cauchy part governs the longitudinal wave propagation. We 
provide explicit examples of the effectiveness of the &4-decomposition. A complete class 
of anisotropic media is proposed that allows pure polarizations in arbitrary directions, 
similarly as in an isotropic medium. 

Key index words: anisotropic elasticity tensor, irreducible decomposition, Cauchy rela- 
tions, null Lagrangians, acoustic tensor 
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1 Introduction and summary of the results 

Consider an arbitrary point P Q , with coordinates x l , in an undeformed body. Then we deform 
the body and the same material point is named P, with coordinates x % . The position of P is 
uniquely determined by its initial position P , that is, x % = x % (x 3 ). The displacement vector 
u = u l di is then defined as 

u l (x{) = x l (O - xl (i, j, ■■• = 1,2,3). (1) 

We distinguish between lower (covariant) and upper (contravariant) indices in order to have 
the freedom to change to arbitrary coordinates if necessary. 

1.1 Strain and stress 

The deformation and the stress state of an elastic body is, within linear elasticity theory, 
described by means of the strain tensor e ^ and the stress tensor a" u . The strain tensor as well 
as the stress tensor are both symmetric, that is, := §(% —Sji) = and cr^ = 0, see Love 
(1927) [25|, Landau & Lifshitz (1986) [22], Haussiihl (2007) [E], Marsden & Hughes (1983) [2B], 
and Podio-Guidugli (2000) [29]; thus, and <7 JJ both have 6 independent components. 
The strain tensor can be expressed in terms of the displacement vector via 

£ij = 9kidju k + g kj diU k =: 2g k{i d j) u k , (2) 

whereupon is the metric of the three-dimensional (3d) Euclidean background space and 
di := d/dx l . In Cartesian coordinates, we have = 2d^Uj), with Wj := gi k u k . 
The stress tensor fulfills the momentum law 

dj(T ij + /;// = pi? ; (3) 

here p is the mass density, b l the body force density, and a dot denotes the time derivative. 

1.2 Elasticity tensor 

The constitutive law in linear elasticity for a homogeneous anisotropic body, the generalized 
Hooke law, postulates a linear relation between the two second-rank tensor fields, the stress 
a 1 ^ and the strain e k i. 

o' :> r" A " (4) 

The elasticity tensor has the physical dimension of a stress, namely force/area. Hence, 
in the International System of Units (SI), the frame components of C'^ kl are measured in 
pascal P, with P:=N/m 2 . 

In 3d, a generic fourth-order tensor has 81 independent components. It can be viewed as 
a generic 9x9 matrix. Since and <7jj are symmetric, certain symmetry relations hold also 
for the elasticity tensor. Thus, 

e M = C^=0. (5) 
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This is the so-called right minor symmetry. Similarly and independently, the symmetry of 
the stress tensor yields the so-called left minor symmetry, 

o-fal = C [ij]kl = . (6) 

Both minor symmetries, (~~~) and ([6]), are assumed to hold simultaneously. Accordingly, the 
tensor C^ kl can be represented by a 6 x 6 matrix with 36 independent components. 

The energy density of a deformed material is expressed as W = § o^Sij. When the Hooke 
law is substituted, this expression takes the form 

W = \ C*%e H . (7) 

The right-hand side of (J7J) involves only those combinations of the elasticity tensor com- 
ponents which are symmetric under permutation of the first and the last pairs of indices. 
In order to prevent the corresponding redundancy in the components of C l ^ kl , the so-called 
major symmetry, 

C ijkl _ C klij = Q (g) 

is assumed. Therefore, the 6x6 matrix becomes symmetric and only 21 independent com- 
ponents of C^ kl are left over. 

As a side remark we mention that the components of a tensor are always measured with 
respect to a local frame e a = e t a di] here a = 1,2,3 numbers the three linearly independent 
legs of this frame, the triad. Dual to this frame is the coframe -& 13 = efd%i, see Schouten 
(1989) [3*5] and Post (1997) [50]. For the elasticity tensor, the components with respect to 
such a local coframe are C a/3lS := ti a ef e^ 1 e\ C l ^ kl . They are called the physical components 
of C. For simplicity, we will not set up a frame formalism since it doesn't provide additional 
insight in the decomposition of the elasticity tensor, that is, we will use coordinate frames <9j 
in the rest of our article. 

Because of (03), (~~~), and (jHJ), we have the following 

Definition: A J^th rank tensor of type (q) qualifies to describe anisotropic elasticity if 

(i) its physical components carry the dimension of force/area (in SI pascal), 

(ii) it obeys the left and right minor symmetries, 

(iii) and it obeys the major symmetry. 

It is then called elasticity tensor (or elasticity or stiffness) and, in general, denoted by . 

In Sees. 2.1 and 2.2 we translate our notation into that of Voigt, see Voigt (1928) [59] and 
Love (1927) [215], and discuss the corresponding 21-dimensional vector space of all elasticity 
tensors, see also Del Piero (1979) [T2]. 

Incidentally, in linear electrodynamics, see Post (1962) [50], Hehl & Obukhov (2003) [18] . 
and Itin (2009) [20], we have a 4-dimensional constitutive tensor ^ uaK = —^ u ^ aK = — -^lvko _ 
X aKflu , with fi, v, ... = 0, 1, 2, 3. Surprisingly, this tensor corresponds also to a 6 x 6 matrix, 
like C l i kl in elasticity. The major symmetry is the same, the minor symmetries are those of 
an antisymmetric pair of indices. 
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1.3 Decompositions of the 21 components elasticity tensor 

In Sec. 2.3, we turn first to an algebraic decomposition of C l ^ kl that is frequently discussed in 
the literature: the elasticity tensor is decomposed into the sum of the two tensors M l ^ hl := 

C iUk)i and N m ._ C imi^ which 

are symmetric or antisymmetric in the middle pair of indices, 
respectively. We show that M and N fulfill the major symmetry but not the minor symmetries 
and that they can be further decomposed. Accordingly, this reducible decomposition does 
not correspond to a direct sum decomposition of the vector space defined by C, as we will 
show in detail. Furthermore we show that the vector space of M is 21-dimensional and that 
of N 6-dimensional. 

Often, in calculation within linear elasticity, the tensors M and N emerge. They are 
auxiliary quantities, but due to the lack of the minor symmetries, they are not elasticities. 
Consequently, they cannot be used to characterize a certain material elastically. These quan- 
tities are placeholders that are not suitable for a direct physical interpretation. 

Subsequently, in Sec. 2.4, we study the behavior of the physical components of C under 
the action of the general linear 3d real group GL(3, R). The GL(3,M) commutes with the 
permutations of tensor indices. This fact yields the well-known relation between the action 
of GL(3,M) and the action of the symmetry group S p . Without restricting the generality of 
our considerations, we will choose local coordinate frames di for our considerations. 




Figure 1. A tensor T % i kl of rank 4 in 3- dimensional (3d) space has 3 4 = 81 independent 
components. The 3 dimensions of our image represent this 81 d space. The plane C depicts 
the 21 dimensional subspace of all possible elasticity (or stiffness) tensors. This space is span 
by its irreducible pieces, the 15d space of the totally symmetric elasticity S (a straight line) 
and the 6d space of the difference A = C — S ( also depicted, I am sad to say, as a straight 
line). — Oblique to C is the 21 d space M of the reducible M -tensor and the 6d space of the 
reducible N -tensor. The C "plane" is the only place where elasticities (stiffnesses) are at 
home. The spaces M and N represent only elasticities, provided the Cauchy relations are 
fulfilled. Then, A = N = and M and C cut in the 15d space of S. Notice that the spaces 
M and C are intersecting exactly on S. 
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In this way, we arrive at an alternative decomposition of C into two pieces S and A, which 
is irreducible under the action of the GL(3, R). The known device to study the action of S p 
is provided by Young's tableaux technique. For the sake of completeness, we will be briefly 
describe Young's technique in the Appendix. In an earlier paper, see Hehl & Itin (2002) [17], 
we discussed this problem already, but here we will present rigorous proofs of all aspects of 
this irreducible decomposition. It turns our that the space of the <S-tensor is 15-dimensional 
and that of the A-tensor 6- dimensional. 

In Sec. 2.5, we compare the reducible MTV-decomposition of Sec. 2.3 with the irreducible 
SA-decomposition of Sec. 2.4 and will show that the latter one is definitely to be preferred 
from a physical point of view. The formulas for the transition between the MN- and the SA- 
decomposition are collected in the Propositions 9 and 11. The irreducible 4th rank tensor A^ kl 
can alternatively be represented by a symmetric 2nd rank tensor Ay (Proposition 10). We 
visualized our main results with respect to the reducible and the irreducible decompositions 
of the elasticity tensor in Figure 1, for details, please see Sec. 2. 5. 

The group GL(3, R), which we are using here, provides the basic, somewhat coarse-grained 
decomposition of the elasticity tensor. For finer types of irreducible decompositions under 
the orthogonal subgroup SO (3) of GL(3,R), see, for example, Rychlewski (1984) [21], Walpole 
(1984) @U, Surrel (1993) [351, Xiao (1998) [33], and the related work of Backus (1970)0 and 
Baerheim (1993) [3]. As a result of our presentation, the great number of invariants, which 
emerge in the latter case, can be organized into two subsets: one related to the S piece and 
the other one to the A piece if C. 

1.4 Physical applications and examples 

Having now the irreducible decomposition of C at our command, Sec. 3 will be devoted to 
the physical applications. In Sec. 3.1, we discuss the Cauchy relations and show that they 
correspond to the vanishing of one irreducible piece, namely to A = or, equivalently, to 
A = 0. As a consequence, the totally symmetric piece S can be called the Cauchy part of 
the elasticity tensor C, whereas the A piece measures the deviation from this Cauchy part; 
it is the non-Cauchy part of C. The reducible pieces M and N defy such an interpretation 
and are not useful for applications in physics. In this sense, we can speak of two kinds of 
elasticity, a Cauchy type and a non-Cauchy type. 

In Sec. 3.2, this picture is brought to the elastic energy and the latter decomposed in a 
Cauchy part and its excess, the "non-Cauchy part" . In other words, this distinction between 
two kinds of elasticity is also reflected in the properties of the elastic energy. 

A null Lagrangians is such an Lagrangian whose Euler-Lagrange expression vanishes, see 
Crampin & Saunders (2005) [IT]; we also speak of a "pure divergence". In Sec. 3.3, elastic null 
Lagrangians are addressed, and we critically evaluate the literature. We show — in contrast 
to a seemingly widely held view — that for an arbitrary anisotropic medium there doesn't 
exist an elastic null Lagrangian. The expressions offered in the corresponding literature are 
worthless as null Lagrangians, since they still depend on all components of the elasticity 
tensor. We collect these results in Proposition 12. 

In Sec. 3.4, we define for acoustic wave propagation the Cauchy and non-Cauchy parts 
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of the Christoffel (or acoustic) tensor P- 7 = C kl ^n^nij p (n^ = unit wave covector, p = 
mass density). We find some interesting and novel results for the Christoffel tensor(, see 
Propositions 13 and 14). In Sec. 3.5, we investigate the polarizations of the elastic wave. We 
show that the longitudinal wave propagation is completely determined by the Cauchy part 
of the Christoffel tensor, see Proposition 15. In Proposition 16 a new result is presented on 
the propagation of purely polarized waves; we were led to these investigations by following 
up some ideas about the interrelationship of the symmetry of the elasticity tensor and the 
Christoffel tensor in the papers of Alshits and Lothe (2004) pQ and Bona et al. (2004, 2007, 
2010)0 El C]. 

In Sec. 4, we investigate examples, namely isotropic media (Sec. 4.1) and media with 
cubic symmetry (Sec. 4.2). Modern technology allows modeling composite materials with 
their effective elastic properties, see Tadmore & Miller (2011) [37]. Irreducible decomposition 
of the elasticity tensor can be used as a guiding framework for prediction of certain features 
of these materials. As an example, we presented in Proposition 17 a complete new class 
of anisotropic materials that allow pure polarizations to propagate in arbitrary directions, 
similarly as in isotropic materials. 

1.5 Notation 

We use here tensor analysis in 3d Euclidean space with explicit index notation, see Sokolnikoff 
(1951) [31] and Schouten (1954, 1989) [321 133]- Coordinate (holonomic) indices are denoted by 
Latin letters i, j, k, . . . ; they run over 1, 2, 3. Since we allow arbitrary curvilinear coordinates, 
covariant and contravariant indices are used, that is, those in lower and in upper position, 
respectively, see Schouten (1989) [33]. Summation over repeated indices is understood. We 
abbreviate symmetrization and antisymmetrization over p indices as follows: 

(k H ■■■ i p ) 
[k h ■■■ ip] 

The Levi-Civita symbol is given by = +1, —1,0, for even, odd, and no permutation of 
the indices 123, respectively; the analogous is valid for e y . The metric ds 2 = gijdx 1 ® dx^ 
has Euclidean signature. We can raise and lower indices with the help of the metric. In 
linear elasticity theory, tensor analysis is used, for example, in Love (1927) [25], Sokolnikoff 
(1956) [35], Landau & Lifshitz (1986) [22], and Haussiihl (2007) [16]. For a modern presenta- 
tion of tensors as linear maps between corresponding vector spaces, see Marsden & Hughes 
(1983) [26], Podio-Guidugli (2000) [29], and Hetnarski & Ignaczak (2011) p]. 

In Hehl & Itin (2002) [17], we denoted the elasticity moduli differently. The quantities 
of our old paper translate as follows into the present one: Wc^ kl = S l ^ kl , ( 2 )c tJ ' w = A^ kl , 

(2)£ijkl = jyijkl^ an( J c i(jk)l ^ M^ kl . 



~|{^^ over & U permutations of i\ %2 ... ip] , 

— over all even perms. — over all odd perms.} . 



Y. Itin and F.W. Hehl 



The constitutive tensor of linear elasticity 



7 



2 Algebra of the decompositions of the elasticity tensor 
2.1 Elasticity tensor in Voigt's notation 

The standard "shorthand" notation of (7 y is due to Voigt, see Voigt (1928) [39] and Love 
(1927) [25]. One identifies a symmetric pair {ij} of 3d indices with a multi-index / that has 
the range from 1 to 6: 

11 -> 1 , 22 -»■ 2 , 33 3 , 23 -»■ 4 , 31 ->• 5 , 12 6 . (9) 

Then the elasticity tensor is expressed as a symmetric 6x6 matrix C IJ . Voigt's notation is 
only applicable since the minor symmetries (JSJ) and ([6]) are valid. Due to the major symmetry, 
this matrix is symmetric, = 0. Explicitly, we have 



£1111 


£1122 


£1133 


£1123 


£1131 


£1112 _ 




-£ll 


£12 


£13 


£14 


£15 


£16 


* 


£2222 


£2233 


£2223 


£2231 


£2212 




* 


£22 


£23 


£24 


£25 


£26 


* 


* 


£3333 


£3323 


£3331 


£3312 




* 


* 


£33 


£34 


£35 


£36 


* 


* 


* 


£2323 


£2331 


£2312 




* 


* 


* 


£44 


£45 


£46 


* 


* 


* 


* 


£3131 


£3112 




* 


* 


* 


* 


£55 


£56 




* 


* 


* 


* 


£1212 




* 


* 


* 


* 


* 


£66 



For general anisotropic materials, all the displayed components are nonzero and independent 
of one another. The stars in both matrices denote those entries that are dependent due to 
the symmetries of the matrices. 



2.2 Vector space of the elasticity tensor 

The set of all generic elasticity tensors, that is, of all fourth rank tensors with the minor 
and major symmetries, builds up a vector space. Indeed, a linear combination of two such 
tensors, taken with arbitrary real coefficients, is again a tensor with the same symmetries. 
We denote this vector space by C. 

Proposition 1. For the space of elasticity tensors, 

dimC = 21. (11) 

Proof. The basis of C can be enumerated by the elements of the matrix C IJ . For instance, 
the element C 11 = C 1111 is related to the basis vector 

E!=di®di®di®di; (12) 

the element C 12 = C 1122 corresponds to the basis vector 

E 2 = - (<9i ® di <8> d 2 ® d 2 + d 2 ® d 2 <g> d x <g> fa) ; (13) 
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the element C = C corresponds to the basis vector 

E 3 = r (di ® 9i <8> <9 3 ® <9 3 + <9 3 <g> <9 3 ® <9i <g> ft) . (14) 

To the element C 14 = C 1123 , we relate the basis vector 

E A = - {d 1 ®d 1 ®d 2 ®d 2 + d 1 ®d l ®d 2 ®d 2 + d 2 ®d- i ®d 1 ®d l + d- i ®d 2 ®d 1 ® d x ) . 

(15) 

In this way, a set of vectors {Ex, • • • , E 2 i} is constructed. Since there are no relations between 
the 21 components C IJ , all these vectors are linearly independent. Moreover every elasticity 
tensor can be expanded as a linear combination of {Ex, ■ ■ ■ , E 2 x}. Thus a basis of of the 
vector space C consists of 21 vectors. □ 

2.3 Reducible decomposition of C^ kl 

2.3.1 Definitions of M and N and their symmetries 

In the literature on elasticity, a special decomposition of C 1 ^ into two tensorial parts is 
frequently used, see, for example, Cowin (1989) [9], Campanella & Tonton (1994) [8], Podio- 
Guidugli (2000) [29], Weiner (2002) (42], and Haussuhl (2007) [16]. It is obtained by sym- 
metrization and antisymmetrization of the elasticity tensor with respect to its two middle 
indices: 

M ijkl ._ Qi{jk)l ^ N ijkl ._ Qi\jk]l ^ with Qijkl = M ijkl + N ijkl _ 

Sometimes the same operations are applied for the second and the fourth indices. Due to 
the symmetries of the elasticity tensor, these two procedures are equivalent to one another. 

We recall that the elasticity tensor fulfills the left and right minor symmetries and the 
major symmetry: 

Q[ij]kl _ q (jv[kl\ = Q . Qijkl _ Qklij _ g Q<jA 

Proposition 2. The major symmetry holds for both tensors M % i kl and N^ kl . 

Proof. We formulate the left-hand side of the major symmetries and substitute the definitions 
given in (II 6ft : 

J^ijkl _ J^fklij _ Qi{jk)l _ Qk(li)j _ _ ^(jijkl _|_ Qikjl _ Qklij _ (jkilj^ _ g 
jyijkl _ jyklij _ (ji[jk]l _ (jk[li]j _ _ ^(jijkl _ Qikjl _ Qklij (jkilj^ _ g 

□ 

Proposition 3. In general, the minor symmetries do not hold for the tensors M^ kl and 
jyijkl^ 
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Proof. We formulate the left minor symmetries for M and N and use again the definitions 
from (TT6l) : 

1 1 1 

j^j-[ij]kl _ _ ^p[ij]kl _|_ Q[i\k\j]l^ _ ±(jk[ij]l _ _^kijl ^Q) 
Ar[ij]kl _ £ (fj[ij]kl _ (j[i\k\j]l\ _ __fjk[ij]l _ _\j^kijl . ^l) 

2 v > 2 2 v ; 

here indices that are excluded from the (anti)symmetrization are enclosed by vertical bars. 
Both expressions don't vanish in general. Moreover, we are immediately led to M^ kl = 

_J\f[ij]kl _^ q 

Using the major symmetry of Proposition 2, we recognize that the right minor symmetries 
M ij[ki] = q and N ij[ki] = o do not hold either, since M ij ^ = and = . □ 

Consequently, the tensors M % i kl and N l ^ kl do not belong to the vector space C and cannot 
be written in Voigt's notation. Thus, these partial tensors M and N themselves cannot serve 
as elasticity tensors for any material. 

2.3.2 Vector spaces of the N- and M -tensors 

We will denote the set of all TV-tensors by J\f. It is a vector space. Indeed, N^ kl is defined 
as a fourth rank tensor which is skew-symmetric in the middle indices and constructed from 
the elasticity tensor. A linear combination of such tensors aN l ^ kl + f3N % ^ kl will be also skew- 
symmetric. Moreover, it can be constructed from the tensor aC % ^ kl + j3C l i kl , which satisfies 
the basic symmetries of the elasticity tensor. 

A simplest way to describe a finite dimensional vector space, is to write-down its basis. In 
the case of a tensor vector space, it is enough to enumerate all the independent components 
of the tensor. 

Proposition 4. For the space of N -tensors, dim M = 6. 

Proof. We can write-down explicitly six components of the "antisymmetric" tensor N^ kl as 
N 1122 = I (C 12 - C m ) , iV 1133 = i (C 13 - C 55 ) , N 1123 = ~ (C 14 - C 56 ) , 

^2233 = I ^23 _ ^44) ? ^2231 = l _ ^25 _ ^46) ? ^1233 = l _ ^36 _ ^45) (22) 

2 2 2 



All other components vanish or differ from the components given in ( 1221) only in sign. Since 
all the components C IJ with I < J are assumed to be independent, then the components of 
N % i kl of ( 122]) are also independent. □ 

The set of all M-tensors is also a vector space, which we denote by Ai. 



Proposition 5. For the space of M-tensors, dimA^ = 21. 
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Proof. The dimension of the vector space Ai can also be calculated by considering the inde- 
pendent components of a generic tensor M %jkl . We find, 



71 /r\ 11 1 
M 


^ 1 


n ^1113 


^-y15 a ,r1 1 12 
C , M 


^-ylfi „ villi 

= (7 , M 


/^y22 
= ^ ) 






= c 26 , 


M 3333 = 


C 33 ii/f 3332 


= C 34 , M 3331 


= c 35 , 




M 2331 


= c 45 , 


M 3221 = 


C 46 , M 3113 


= C 55 , M 3112 


= c 56 , 




M 1122 




+ c m ) , 


M 1133 = I 

2 


(C 13 + C 55 ) , 


M 2223 = 


c 24 


M 1123 


- 5 ^ 


+ c 56 ) , 


M 2233 = I 

2 


(C<23 + C 44) ^ 


M 2332 = 


c 44 


M 2231 


= \ (c 25 


+ c 46 ) , 


M 1233 = 1 

2 


(C<36 + ^45) ^ 


M 1221 = 





(23) 

All these 21 components are linearly independent, thus the dimension of the vector space 
Ai is at least 21. However, since every element of Ai is defined in terms of 21 independent 
elastic constants C IJ , the dimension of Ai cannot be greater than 21. □ 

2.3.3 Algebraic properties of M and N tensors 

Observe some principal features of the tensors M l ^ kl and N v : 

(i) Inconsistency. In general, a certain component of C^ kl , say C 1223 , can be expressed 
in different ways in terms of the components of M yM a N^ kl : 

£<1223 M 1223 _|_ ]Y 122 3 ^2312 j^2132 (2J) ^46 

On the other hand, we have 

(J 1223 = (^2123 (22) ^2123 _|_ jy2123 ^5) 



With 

M 2123 = - (C 46 + C 25 ) and A 2123 = -(C 46 -C 25 )^0, (26) 

we recover the result in ( 1241) , but is was achieved with the help of a non- vanishing component 
of N. 

(ii) Reducibility . Since, in general, the tensor M l i kl is not completely symmetric, a finer 
decomposition is possible, 

M ijki = M (ijki) _|_ K ijki _ ^27) 
Accordingly, C* 1 - 7 can be decomposed into three tensorial pieces: 

Qijkl _ j^(ijkl) _|_ j^ijkl _|_ jyij'fcZ ^8) 

(^mj Vector spaces. The "partial" vector spaces, A-4 and A/", are not subspaces of the 
vector space C and their sum Ai + A/" is not equal to C. 

Thus, the MTV-decomposition is problematic from an algebraic point of view. Our aim is 
to present an alternative irreducible decomposition with better algebraic properties. 
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2.4 Irreducible decomposition of 

2.4.1 Definitions of S and A and their symmetries 

In Eg. (1146[) of the Appendix, we decomposed a fourth rank tensor irreducibly. Let us apply 
it to the elasticity tensor C l ^ kl . Since the dimension of 3d space is less than the rank of 
the tensor, the last diagram in ( 114611 . representing C^ kl \ is identically zero. Also the minor 
symmetries remove some of the diagrams. Dropping the diagrams which are antisymmetric 
in the pairs of the indices (ii, 12) and (k, k), we are eventually left with the decomposition 



1-1 



12 <a> z 3 



14 



1 

I 



k 


k 


k 


k 



+ P 



k 


k 


k 


k 





+ 



k 


k 


k 


k 





7 



k 


k 


k 


u 



(29) 



Let us now apply the major symmetry. It can be viewed as pair of simultaneous permutations 
k -H- k an d k "H- k. The first and the last diagrams in (129]) are invariant under those 
transformations. The two diagrams in the middle change their signs and are thus identically 
zero. Thus, we are left with irreducible parts: 

1 



1-2 



1! 



k 


k 


k 


k 



7 



k 


k 


k 


k 



(30) 



In correspondence with the first table, the first subtensor of C ftJ is derived by complete 
symmetrization of its indices 



gijki = (7 + (i u i 2 ) + i 3 ) + ... + (i 1} { 2} i 3 ) + . . . + (i u { 2) i 3j j 4 ))cy fc ' 
Here the parentheses denote the cycles of permutations. Consequently, 



Cjijkl £t(ijkV) 



J_ ( (jijkl _|_ Qjikl _|_ Qkjil _|_ . . . _|_ (jkijl _|_ . . . _|_ Qlijk 

4! V 



(31) 



(32) 



On the right hand side, we have a sum of 4! = 24 terms with all possible orders of the indices. 
If we take into account the symmetries (jSJ), (EL and (jHJ) of C iJ , we can collect the terms: 



gijkl £j(ijkl) _{£fijkl _|_ (jiklj _|_ (jiljk^ 

3 

According to fl30|) . the second irreducible piece of C^ kl can be defined as 

j^ijkl ^jijkl £j(ijkl) ^jijkl gijkl 



(33) 



(34) 



This result can also be derived by evaluating the last diagram in (1301) . Substitution of (1331) 
into the right-hand side of fl34l) yields 

1 



ki 



/^^-jijkl ^jilkj £jiklj\ 



(35) 



If we totally symmetrize the left- and the right-hand sides of (I34|) . an immediate consequence 
is 

A (tjkl) = . (36) 

If we symmetrize (1351) with respect to the indices jkl, we recognize that its right-hand side 
vanishes. Accordingly, we have the 



Y. Itin and F.W. Hehl 



The constitutive tensor of linear elasticity 



12 



Proposition 6. The tensor A fulfills the additional symmetry 

A Ojki) _ g Qr A ijki + A ikij + A ujk = _ ^ 

2.4.2 Vector spaces of the S- and A-tensors 

We denote the 21-dimensional vector space of C by C. The irreducible decomposition of C 
signifies the reduction of C to the direct sum of its two subspaces, S C C for the tensor S, 
and A C C for the tensor A, 

C = S@A. (38) 

The vector spaces S and A have only zero in their intersection and the decomposition of the 
corresponding tensors is unique. According to Proposition 1, the sum of the dimensions of the 
subspaces is, equal to 21. The two irreducible parts S %3 and A ljkl preserve their symmetries 
under arbitrary linear frame transformations. In particular, they fulfill the minor and major 
symmetries of C %3kl likewise. Accordingly, S and A, or S alone (but not A alone, as we will 
see later) can be elasticity tensors for a suitable material — in contrast to M and N. 
The dimensions of the vector spaces of 5* and A can now be easily determined. 

Proposition 7. For the vector space S of the tensors S l i kl , 

dimS = 15. (39) 

Proof. The number of independent components of a totally symmetric tensor of rank p in 
n dimensions is ( n+ P -1 ) = ( n ~^\J~ p ) or, for dimension 3 and rank 4, (^) = 15, see Schouten 
(1954)|32]. □ 

According to Proposition 1, we have dimC = 21. Because of ( J38l) and ( |39l) . we have 
Proposition 8. For the vector space A of the tensors A 13 , 

dirndl = 6. (40) 

2.4.3 Irreducible parts in Voigt's notation 

In Voigt's 6d notation we have 

c u = s u + A u with C [ij] = S [u] = A [u] = _ ( 41 ) 

The 6x6 matrix S IJ has 15 independent components. We choose the following ones [see 
Voigt (1928), Eq. (36) on p.578], 



s n 


= c n , 


g22 _ q22 


s 33 


= C 33 , 




s 15 


= c 15 , 


Cjil6 (J^ 


S 26 


= C 2 \ 




s 2i 


= c 2 \ 


g34 _ £<34 


S 35 


= C 3 \ 




s 12 




+ 2C m ) , 


S 13 = 




+ 2C 55 ) , 






+ 2C 56 ) , 


S 23 = 


\{c 23 


+ 2C 44 ) , 


s 25 


- 1 ^ 


+ 2C 46 ) , 


S 36 = 


l (c - 


+ 2C 45 ) . 



(42) 
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The 6x6 matrix of A IJ has 6 independent components. We choose the following ones, 

A 12 = - (C 12 - C 66 ) , A 13 = - (C 13 - C 55 ) , A u = - (C 14 - C 56 ) , 
3 3 3 

A 23 = - (C 23 - C 44 ) , A 25 = - (C 25 - C 46 ) , A 36 = - (C 36 - C 45 ) . (43) 
3 3 3 

The decomposition (14 ip can be explicitly presented as 



C ll 


c 12 


C 13 




c 15 








5 13 




5 15 5-16- 




* 


c 22 


c 23 


C 2A 


c 25 


c 26 




* 


s 22 




£24 


£25 5 26 




* 


* 


c 33 


C 3A 


c 35 


c 36 




* 


* 


5 33 




5 35 5 36 




* 


* 


* 




c 45 


c 46 




* 


* 


* 


s 23 


s 36 s 25 




* 


* 


* 


* 




c 56 




* 


* 


* 


* 


s 13 s u 




* 


* 


* 


* 


* 






* 


* 


* 


* 


* s 12 _ 
















"0 


A 12 


A 13 


A 14 




















* 





A 23 





A 25 

















* 


* 











A 36 












+ 


* 


* 


* 


-kA 23 


-IA 36 
















* 


* 


* 


* 


-\A A3 
















* 


* 


* 


* 


* 





(44) 



Here, we use boldface for the independent components of the tensors. Note that all three 
matrices are symmetric. 



2.5 Comparing the SA- and the MTV-decompositions with each 
other 

We would now like to compare the two different decompositions: 

(jijkl _ gijkl _|_ j^ijkl _ j^ijkl _|_ jyijkl 
21 15 6 21 6 

The dimensions of the corresponding vector spaces are displayed explicitly. This makes it im- 
mediately clear that A can be expressed in terms of N and vice versa. Take the antisymmetric 
part of (j45p with respect to j and k and find: 

Proposition 9. The auxiliary quantity can be expressed in terms of the irreducible 

elasticity A l ^ kl as follows: 

jyijkl _ J^i[jk]l _ ^g-j 

Its inverse reads, 

j^ijkl _ 4 j^ij{kl) ^ ^ij-j 
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Proof. Resolve f )46p with respect to A. For this purpose we recall that A obeys the right 
minor symmetry: A % ^ kl = A l ^ kl \ This suggests to take the symmetric part of ()46p with 
respect to k and I. Then, 

jyij{M) _ i ^ijki _ j^,(ki)j^ _ i ^2A^ kl — A ikl i — A ilk ^ 

i ^2A^ kl — (— A il ^ k — A 1 ^ 1 ) — A ilk ^ = iA^ kl (48) 

□ 

Both, A and N have only 6 independent components. In 3d this means that it must be 
possible to represent them as a symmetric tensor of 2nd rank. With the operator |e m ij, we 
can always map an antisymmetric index pair ij to a corresponding vector index m. The 
tensor A l i kl has 4 indices, that is, we have to apply the e operator twice. Since A l ^ kl obeys 
the left and right minor symmetries, that is, A^ kl = A l ^ kVi = 0, the e has always to transvect 
one index of the first pair and one index of the second pair. This leads, apart from trivial 
rearrangements, to a suitable definition. 

Proposition 10. The irreducible elasticity A l i kl can be equivalently described by a symmetric 
2nd rank tensor 

A mn ■ "T^mik^njlA ^ , (^9) 

with the inverse 

N ijlk = e ikm e jln Amn Qr A ijkl = l e i(fc|"V'l' )n A mn . (50) 

3 

Proof. The symmetry of A mn can be readily established: 

^[m.n] = g (tmiktnjlA J — € n ik€ m jlA^ ) = -^raik^njli^A^ — A^ ) = 0. (51) 

Eq. ( l50l) i can be derived by substituting ( 1491) into its right-hand side and taking care of ( j46l ). 
Eq. ( 1501) '? then follows by applying ( I4"T1) . □ 

The symmetric 2nd rank tensor A mn (differing by a factor 2) was introduced by Haussiihl 
(1983, 2007), its relation to the irreducible piece A was found by Hehl & Itin (2002). By 
means of fl34l) . A mn can be calculated directly from the undecomposed elasticity tensor: 

A mn = —Grnik^njl C ' J . (52) 

The tensor M, like N, can also be expressed in terms of irreducible pieces: We substitute 
(gSD into g5]), 

Qijkl _ gijkl _|_ j^ijkl _ J^ijkl _|_ j^i[jk]l ^ggj 

and resolve it with respect to M: 

J^ijkl _ gijkl _|_ j^ijkl _ j^i[jk]l 

= S ljkl + i {2A i3kl - A ljkl + A ikjl ) . (54) 

We collect the terms and find 
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Proposition 11. The auxiliary quantity M y can be expressed in terms of the irreducible 
elasticities as follows: 

M ijki _ s m + A i(jk)i _ ^ 

According to Proposition 3, the reducible parts M and N don't obey the left and the 
right minor symmetries. Consequently — in contrast to S and A, which both obey the major 
and the minor symmetries — M and N cannot be interpreted as directly observable elasticity 
tensors. Therefore, in all physical applications we are forced to use eventually the irreducible 
pieces S and A. The reducible parts M and N are no full-fledged substitutes for them and 
can at most been used for book keeping. 



3 Physical applications of the irreducible decomposi- 
tion 

3.1 Cauchy relations, two kinds of elasticity 

Having the SA-decomposition at our disposal, it is clear that we can now classify elastic 
materials. The anisotropic material with the highest symmetry is that for which the A 
elasticity vanishes: 

A ijkl = or A mn = or N ijkl = . (56) 
The last equation, in accordance with the definition (TIB"]) of N, can also be rewritten as 

Qijkl _ Qikjl 

These are the so-called Cauchy relations, for their history, see Todhunter (1960) [58]. The rep- 
resentation in (1571) is widely used in elasticity literature, see, for example, Haussuhl (1983) [15"]. 
Cowin (1989) [9], Cowin & Mehrabadi (1992) [10], Campanella & Tonon (1994)0. Podio- 
Guidugli (2000) [29], Weiner (2002) [12], and Hehl & Itin (2002) pi]. 

In Voigt's notation, we can use A 13 = in fj4*3|) and find the following form of the Cauchy 
relations, 

C 23 = C 4 \ C 25 = C 4 \ C 36 = C 45 , (58) 

see Love (1927) [25] and Voigt (1928) [39]. Of course, the same result can also be read off from 
d22D for iV = 0. 

Let us first notice that for most materials the Cauchy relations do not hold even approx- 
imately. In fact, the elasticity of a generic anisotropic material is described by the whole set 
of the 21 independent components C IJ and not by a restricted set of 15 independent com- 
ponents obeying the Cauchy relations, see Haussuhl (2007) [IB]. This fact seems to nullify 
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the importance of the Cauchy relations for modern solid state theory and leave them only as 
historical artifact. 

However, a lattice-theoretical approach to the elastic constants shows, see Leibfried 
(1962) [23], that the Cauchy relations are valid provided (i) the interaction forces between 
the atoms or molecules of a crystal are central forces, as, for instance, in rock salt, (ii) each 
atom or molecule is a center of symmetry, and (iii) the interaction forces between the build- 
ing blocks of a crystal can be well approximated by a harmonic potential, see also Perrin 
(1979) [28]. In most elastic bodies this is not fulfilled at all, see the detailed discussion in 
Haussiihl (2007) [16]. Accordingly, a study of the violations of the Cauchy relations yields 
important information about the intermolecular forces of elastic bodies. One should look for 
the deviation of the elasticity tensor from its Cauchy part. Recently, Elcoro & Etxebarria 
(2011) [13] pointed out that the situation is more complex than thought previously, for details 
we refer to their article. 

This deviation measure, being a macroscopic characteristic of the material, delivers im- 
portant information about the microscopic structure of the material. It must be defined in 
terms of a unique proper decomposition of the elasticity tensor. Apparently, the tensor N l i kl 
cannot serve as a such deviation, in contrast to the stipulations of Podio-Guidugli (2000), 
for example, because it is not an elasticity tensor and because its co-partner M !jH has 21 
components, that is, as many as the elasticity tensor C 1 ^ itself. Only when N l ^ kl = is 
assumed, the tensor M^ kl is restricted to 15 independent components. The problem of the 
identification of the deviation part is solved when the irreducible decomposition is used. 
In this case, we can define the main or Cauchy part, as given by the tensor S l ^ kl with 15 
independent components, and the deviation or non-Cauchy part, presented by the tensor 
A^ kl with 6 independent components. 

We can view the elasticity tensor as being composed of two independent parts, 5* and 
A. Due to the irreducible decomposition the set of elastic constants can be separated into 
two subsets which are components of two independent tensors, C = S © A. In this way, the 
completely symmetric Cauchy part of the elasticity tensor S l ^ kl has an independent meaning. 
The additional part A^ kl can be referred than as as a non-Cauchy part. Why do pure 
Cauchy materials (S ^ 0, A = 0) and pure non-Cauchy materials (S = 0, A ^ 0) not exist 
in nature? Can such pure types of materials be designed artificially? Or does some principal 
fact forbid the existence of pure Cauchy and pure non-Cauchy materials? These questions 
seem to be important for elasticity theory and even for modern material technology. Note 
that in the framework of the reducible M/V-decomposition such questions cannot even be 
raised. It is because the reducible M and N parts themselves cannot serve as independent 
elasticity tensors. We will address subsequently the reason why pure non-Cauchy materials 
are forbidden. 

3.2 Elastic energy 

In linear elasticity, using the generalized Hooke law, the elastic energy is given by 



(59) 
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Because of the irreducible decomposition C = S + A, we can split this energy in a Cauchy 
part and a non-Cauchy part: 

W = {C W + {nC) W , with '% := -S ijkl e i:j e kl , (nC) W := -A^e^u . (60) 

As we have seen in the Sec. 3.1, this splitting makes perfectly good sense in physics, since S 
as well as A, due to their symmetries and dimensions, are elasticity tensors themselves. 

Since the strain be expressed in terms of the displacement gradients according to 

Eij = 2<9(jWj), we can do the analogous for the elastic energy: 

W = ^C ijkl £ij E kl = 2C tjkl d {lUj) d {kUl) = 2C l ' kl d lUj d kUl . (61) 

We can drop both pairs of parentheses q because the corresponding symmetries are imprinted 
already in the elasticity tensor. 

3.3 Null Lagrangians in linear elasticity? 

We would like now to discuss a proposal on null Lagrangians by Podio-Guidugli (2000) [29]. 
see also literature cited by him. We substitute the decomposition ( fl6l) into ( |6"T1) : 

W = 2C ijkl d iUj d k ui = 2M ijkl d i u j d k u l + 2N ijkl d i u j d k u l . (62) 

We turn our attention now to the last term and to the antisymmetry of N, namely to 
N i(jk)i = ^ There 

is a subtlety involved. We want to integrate partially in the last term. 
For exploiting the antisymmetry of iV in j and k, we need the partial differentials dj and d k . 
Therefore, we use the left minor symmetry of C and rewrite (!62|) as 

W = 2C ijkl d jUi d k ui = 2M ljkl d j u i d k u l + 2N ijkl d jUl d k ui . (63) 

Of course, the sum of the two terms has not changed, but each single term did change since 
M and TV do not obey the left minor symmetry. Now we can partially integrate in the last 
term: 

2N ljkl d jUi d kUl = 2N ijkl [dj { Ui d k ui) - Uidjdkm] . (64) 
Since N 1 ^^ 1 = and dyd k ] = 0, the last term drops out and we are left with 

W = 2M ijkl d jUl d kUl + 2N ijkl d j {u l d k u l ) . (65) 

The second iV-term is a total derivative, as was already shown by Lancia et al.(1995)|21j. 
Thus, only the M-term is involved in the variational principle for deriving the equations of 
motion. This result was interpreted by Lancia et al. (1995) [21] as solving the null Lagrangian 
problem for linear elasticity theory, since the additional M-tensor is the only one that is 
involved in the equilibrium equation. However, this statement does not have an invariant 
meaning. Indeed, the remaining part with the M-tensor has exactly the same set of 21 
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independent components as the initial C-tensor, see Proposition 5. Moreover, this tensor 
allows a successive decomposition which can include additional null Lagrangian terms. 
We find from ( 165]) . using the Propositions 9 and 11 and subsequently Proposition 6, 

W = {2S llkl - A iljk ) djUidkm + 2A l ^ l d j (uidkUi) . (66) 

Consequently the A-tensor is included in the total derivative term, which vanishes if the 
Cauchy relations hold. However the same A-tensor appears together with the S'-tensor also 
in the first part of the energy functional. 

For static configurations, (159 p can play a role of a Lagrangian functional whose varia- 
tion with respect to the displacement field generates the equilibrium equation. The null- 
Lagrangian is defined as that part of the strain energy functional that does not contribute to 
the equilibrium equation. The problem is to identify the null-Lagrangian part of the elasticity 
Lagrangian and consequently to establish which set of the independent components of the 
elasticity tensor contributes to the equilibrium equation. 

As a cross-check for our considerations, we determine the equilibrium conditions for the 
Lagrangian (16"6"|) . Up to a total derivative term, the variation of the Lagrangian (16"6"|) reads 

SW = (2S ljkl - A lljk ) {djUidkSui + djS Ul d k ui) . (67) 

Since the minor and the major symmetries hold for the S and A tensors, the last two terms 
can be summed up: 

5W = 2 (2S ijkl - A tljk ) djSuidkUi . (68) 

We integrate partially, 

SW = 2 (2S ljM - A iljk ) dj (Suid k ui) - 2 (2S ijkl - A ll]k ) Sutd^m , (69) 

and can read off the equilibrium conditions as 

(2S i]kl - A iljk ) djdkUt = . (70) 

Only at a first glance, this equation seems to be new. Indeed, we can rewrite it by using 
the Propositions 6 and 11, 

(2S 4jU + A i3k[ + A iklj ) dfitui = 2 (S ijkl + A l(lk)l ) djd k ui = 2M ijkl djd kUl = , (71) 

Since M % i kl = C l ^ k ^ 1 , we can substitute it and, because of dydk] = 0, we have the standard 
equilibrium equation 

C^ kl d,d kUl = 0. (72) 
Our calculations confirms the following: 

Proposition 12. Any total derivative term in the elastic energy functional can be regarded 
as only a formal expression. It does not remove any subset of the elastic constants from the 
equilibrium equation, that is, for an arbitrary anisotropic material an elastic null Lagrangian 
does not exist. 
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3.4 Wave equation 

The wave propagation in linear elasticity for anisotropic media is described by the following 
equation: 

pg a u l -C llkl d j d k u l = 0. (73) 

Here the displacement covector u\ is assumed to be a function of the time coordinate and 
of the position of a point in the medium. All other coefficients, the mass density p, the 
elasticity tensor C l ^ kl , and the metric tensor g d , are assumed to be constant; moreover, we 
use Cartesian coordinates, that is, the Euclidean metric reads = diag(l, 1,1). 
We make a plane wave ansatz, with the notation as in Nayfeh (1985) [27]: 



Ui 



Here Ui is the covector of a complex constant amplitude, ( the wave-number, rij the propa- 
gation unit covector, u the angular frequency, and i 2 = — 1. Substituting (|74p into (175]) . we 
obtain a system of three homogeneous algebraic equations 

(puj 2 g u -C i i kl C 2 n j n k )U l = [). (75) 

It has a non-trivial solution if and only if the characteristic equation holds, 

det (pu 2 g il - C ijkl C 2 n jnk ) = . (76) 

With the definitions of the Christoffel tensor 

r u ._ I c<;w (77) 
P 

and of the phase velocity v := oo/C, the system ( 1751) takes the form 

(v 2 g il - r<) Ui = Q, (78) 

whereas the characteristic equation reads 

det (v 2 g il - T il ) = . (79) 

Due to the minor and major symmetries of the elasticity tensor, the Christoffel tensor turns 
out to be symmetric 

T [ij] = . (80) 

For a symmetric matrix, the characteristic equation has only real solutions. Every positive 
real solution corresponds to an acoustic wave propagating in the direction of the wave covector 
Hi. Thus, in general, for a given propagation covector n;, three different acoustic waves are 
possible. The cases with zero solutions (null modes) and negative solutions (standing waves) 
must be treated as unphysical because they do not satisfy the causality requirements. 

A symmetric tensor, by itself, cannot be decomposed directly under the action of the 
group GL(3,R). However, the irreducible decomposition of the elasticity tensor generates 
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the corresponding decomposition of the Christoffel tensor. Substituting the irreducible SA- 
decomposition into (J77J) , we obtain 

V il = S il + A il , (81) 

where 

s u ._ S ijki n . nk= s n ^ and A u = A m n . nk= A n . ( 82 ) 

These two symmetric tensors correspond to the Cauchy and non-Cauchy parts of the elasticity 
tensor. We will call «S y Cauchy Christoffel tensor and A^ non-Cauchy Christoffel tensor. 
Substituting (|33|) and (I35|) into (1521) . we find the explicit expressions 

S il = —(C ijkl + C mj + C Ujk ) nj n k (83) 
3p 

and 

A il = — (2C ljkl - C lMl - C ilkj ) rijrik . (84) 
3p 

If we transvect ni with A d , we observe a generic fact: 

Proposition 13. For each elasticity tensor C^ kl and for each wave covector n^, 

A il n t = . (85) 

(37) 

Proof. We have A ll n\ = A^n^n^ni = A l ^ kV> njUkni = 0. □ 

Proposition 14. The determinant of the non-Cauchy Christoffel tensor is equal to zero, 

det(^')=0. (86) 

Proof. Equation (J85l) can be viewed as a linear relation between the rows of the matrix A d ; 
this proves its singularity. □ 

The wave equation (1781 can be rewritten as 

(v 2 g tl - S il - A il ) Ui = , (87) 

with the characteristic equation 

det (v 2 g il - S il - A 11 ) = . (88) 

We can now recognize the reason why a pure non-Cauchy medium is forbidden. In the case 
S d = 0, the characteristic equation ( )88l) takes the form 



det (v 2 g il - A 11 ) = . (89) 

Since det(A d ) = 0, at least one of its eigenvalues is zero. Such a null- mode wave is forbidden 
because of causality reasons. 

The Christoffel tensor is real and symmetric, thus all its eigenvalues are real and the 
associated eigenvectors are orthogonal. In order to have three real positive eigenvalues, we 
need to satisfy the condition of positive definiteness of the matrix P J . In this case, the 
following possibilities arise: 
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(i) All the eigenvalues are distinct v\ > v\ > v 2 , 

(ii) two eigenvalues coincide v\ > v \ = v 2 or vf — v 2 > v 2 , or 

(iii) three eigenvalues coincide v\ = v 2 = v 2 . 

The Christoffel matrix depends on the propagation vector r u (n), thus the conditions 
indicated above determine the directions in which the wave can propagate. 

3.5 Polarizations 

Equation flTHj) . or its decomposed form flHTj) . represents acoustic wave propagation in an 
elastic medium. It is an eigenvector problem in which eigenvalues v 2 are the solutions of 
fl88l) . In general, three distinct real positive solutions correspond to three independent waves 
Wf/j, ( 2 )[/;, and ^Ui, called acoustic polarizations, see Nayfeh (1985)[2"Tj. 

For isotropic materials, there are three pure polarizations: one longitudinal or compression 
wave with 

U x n = , (90) 

that is, the polarization is directed along the propagation vector, and two transverse or shear 
waves with 

U-n = 0, (91) 

that is, the polarization is normal to the the propagation vector. In general, for anisotropic 
materials, three pure modes do not exist. The identification of the pure modes and the 
condition for their existence is an interesting problem. 

Let us see how the irreducible decomposition of the elasticity tensor, which we applied 
to the Christoffel tensor, can be used here. For a chosen direction vector n, we introduce a 
vector and a scalar according to 

S i := V ij rij , S := T i3 n inj . (92) 

This notation is consistent since S l and S, due to ( 185|) . depend only on the Cauchy part of 
the elasticity tensor, 

S l = S ij nj , S = S ij n inj . (93) 

Proposition 15. Let n % denotes an allowed direction for the propagation of a compression 
wave. Then the velocity of this wave in the direction of n l is determined only by the 
Cauchy part of the elasticity tensor: 

v L = y/S. (94) 
Proof. For the longitudinal wave, Uj = anj. Thus, ( !7H|) becomes 

(v 2 g ij - nj = or v 2 g l3 nj = S ij nj . (95) 



Transvecting both sides of the last equation with rij, we obtain (1941) . 



□ 
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Let us now discuss in which directions the three pure polarizations can propagate. 

Proposition 16. For a medium with a given elasticity tensor, all three purely polarized waves 
(one longitudinal and two transverse) can propagate in the direction n if and only if 

S i = Sn l . (96) 

Proof. Since the Christoffel matrix is symmetric, it has real eigenvalues and three orthogonal 
eigenvectors. We have three pure polarizations if and only if one of these eigenvectors points 
in the direction of n. Consequently, 

T ij nj = vln 1 . (97) 
Since T^nj = S l ^nj = S\ we have S l = v\n % or, after substituting (1911) . we obtain (I9"6"i) . □ 

Accordingly, for a given medium, the directions of the purely polarized waves depend on 
the Cauchy part of the elasticity tensor alone. In other words, two materials, with the same 
Cauchy parts S %3kl of the elasticity tensor but different non-Cauchy parts A l ^ kl , have the same 
pure wave propagation directions and the same longitudinal velocity. 



4 Examples 

4.1 Isotropic media 

In order to clarify the results discussed above, consider as an example an isotropic elastic 
medium. Then, the elasticity tensor can be expressed in terms of the metric tensor g l i as 

with the Lame moduli A and see Marsden & Hughes (1983) [26]. The irreducible decom- 
position of this expression involves the terms 

S ijkl = (A + 2/i) = (<7 V + 9 lk 9 l] + 9 ll 9 ]k ) , (99) 

and 

A ijkl = (2g lj g kl - g lk g lj - g il g jk ) . (100) 

The symmetric second rank tensor A of equation fH91 . which is equivalent to A, reads 

A mn = \e mik e njl A^ kl = ^^e mik e njl (2g^g kl - g ik g 1 ^ - g a g> k ) (101) 

or 

A lJ = ^V (102) 
Of course, (I102f) is much more compact than (11001) . 
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Accordingly, the irreducible decomposition reveals two fundamental combinations of the 
isotropic elasticity constants. We denote them by 

Then we can rewrite the canonical representation (158]) of the elasticity tensor as 

C*'« = a {g^g kl + g* k g l i + g u gi k ) + ft (2<? V - V ~ <?V' fc ) 

= {a + 2(3) g ij g kl + 2 (a - 0) g i{k g l)j . (104) 

The alternative M/V-decomposition can also be derived straightforwardly. We can read 
off from 098p directly that 

M ijki = C i(jk)i = ( A + ^)g^ g k ) 1 + . (105) 

In a similar way we can find for N^ kl the expression 

N ijki = C iim = (A - ^) g^ g k V . (106) 

The key fact is here that this decomposition is characterized by three parameters 

«=^— ? P =/i, ^ = ^~ • ( 107 ) 

In both decompositions, the Cauchy relations are described by the equation 

A = /i. (108) 

Only in the special case when the Cauchy relations hold, the tensors S l ^ kl and M^ kl are 
equal. 

Consider now the energy functional and the problem of the identification of the null 
Lagrangian. When the SA-decomposition is substituted into the energy functional, we remain 
with two irreducible parts with their leading coefficients a and /3. Thus, the number of 
parameters is not reduced. 

As for the M/V-decomposition, the iV term, with its leading coefficient 7', being a total 
derivative, does not contribute to the equilibrium equation. We still remain with two terms 
with their leading parameters a' and 0'. Moreover, as it was shown above, these terms can be 
recovered in the initial functional. It proves that there is no such thing as a null Lagrangian 
in elasticity, not even in the simplest isotropic case. 

In this isotropic case, the main difference between the SA- and the M/V-decomposition be- 
comes manifest. The irreducible S!A-decomposition dictates the existence of two independent 
fundamental parameters of the isotropic medium, whereas the reducible M/V-decomposition 
deals with three linearly dependent parameters. 

Let us calculate now the characteristic velocities of the acoustic waves. Two two Christof- 
fel matrices are 

= a (g i] + 2nW) . (109) 
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and 



Hence, 



A 



f3 (g ij - nW) 



S l = 3an l 



S = 3a 



;no) 



111^ 



Observe that the condition f PB|) of pure polarization is satisfied now identically. Thus, we 
recover the well known fact that in isotropic media every direction allows propagation of 
purely polarized waves. 

In terms of the moduli a and /3, the characteristic equation for the acoustic waves takes 
the form 

det [{y 2 -a + 0)g ij - (2a + /3)n l n j ] = . (112) 
The longitudinal wave velocity 



S = 3a = A — 2ji , 



(113) 



is one solution for this equation. Indeed, in this case, the characteristic equation (|112p turns 



into the identity det(g 



ij 



n n 



0, which proves that (j!13p is an eigenvalue. 



We immediately find additional solutions for this equation. For 



a 



/3 — fj, 



(114) 



f)112p turns into the identity det(n*n J ) = 0. Moreover, the rank of this matrix is equal to 
one: rank(n*n J ) = 1. Consequently, the multiplicity of the root (I114p is equal to two and the 
third eigenvalue is the same one: 



a 



P = V 



115) 



With the use of the irreducible decomposition, we can answer now the following question: 
Does there exist an anisotropic medium in which purely polarized waves can propagate in 
every direction? It is clear that media with the same Cauchy part of the elasticity tensor 
(and arbitrary non-Cauchy part) have the same property in this respect. As a consequence 
we have 

Proposition 17. The most general type of an anisotropic medium that allows propagation 
of purely polarized waves in an arbitrary direction has an elasticity tensor of the form 

a a/3 + 2 Pl a/3 + 2p 2 2p 3 

* a a/3 + 2p 4 2p 5 

* * a 2p 6 

* * * a/3- pi -p 6 —p5 

* * * * a/3 — P 2 — P 3 

* * * * * a/3 — p% 



C ij 



kl 



;ii6) 



where pi, • • • , Pe are arbitrary parameters, 
medium is t>L = y/3~a = y/X — 2p. 



The velocity of the longitudinal waves in this 
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4.2 Cubic media 



Cubic crystals are described by three independent elasticity constants. In a properly chosen 
coordinate system, they can be put, see Nayfeh (1985) [27], into the following Voigt matrix: 



(117) 



£1111 


£1122 


£1133 


£1123 


£1131 


£1112" 




£11 


£12 


£12 
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£2222 


£2233 


£2223 


£2231 


£2212 
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£11 


£12 
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£3333 


£3323 


£3331 


£3312 




* 


* 


£11 
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£2323 


£2331 


£2312 
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£66 
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£3131 


£3112 
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£66 
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£1212 
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* 


£66 



We decompose it irreducibly by using ( H2"j) and ( H3"j) and find the Cauchy part 
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non- Cauchy part 
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c 



11 



p 



(C 12 + 2C 66 ) 
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7 



(£12 _ £66) 



(119) 



Accordingly, the elasticity tensor is expressed in terms of the three new elastic constants a, 
(3, and 7. For the Cauchy relation we have 

£12 _ £66 



7 = 0, or 

The Cauchy part of the Christoffel tenor takes the form 



S' 



an\ + (3(nl + n\) 2j3n 1 n 2 

2f3niri2 an?, + f3{n\ + n 2 ) 

2f3n l n 3 2/3n 2 n 3 



2/3nin 3 
2/3n 2 n 3 
an\ + P{n\ + r?|)_ 



(120) 



'121) 



The corresponding vector and scalar invariants fl93l) read 



S' 



m(an{ + 3/3(?tJ + n§)) 
n 2 (an 2 + 3(3 (nj + n 2 )) 
n 3 {anl + 3f3(n\ + n 2 2 ))_ 



n x ({a-3p)n\ + 3/3) 
n 2 ((d - 3(3)n 2 + 3/3) 
n 3 ((S-30)n§ + 30). 



:i22i 
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and 

S = (a- 3/3)K + n t + n t) + 3/3 , 
respectively. The non-Cauchy part of the Christoffel tensor turns out to be 



.4' 



7 



""2 

n x n 2 
n x n 3 



n x n 2 

9 

n 2 n 3 



nin 3 
n 2 n 3 

2 

-ni 



(123) 



(124) 



The corresponding invariants A^rij and A^riiUj are zero. 

Let us use the pure polarization condition S" 1 = Sn l in order to derive the purely polarized 
propagation directions. We can easily derive all solutions of this equation 



Edges: n\ = 1 , n 2 = n 3 = , etc. The longitudinal velocity is 

vl = 



Face diagonals: n\ — n 2 — \\/2 ,n 3 = ,etc. The longitudinal velocity 



is 



-J2{a + 3f3) = -^2(C ll + C l2 + 2C™). 



(125) 



(126) 



Space diagonals: ri\=n 2 = n 3 = \\^3,etc. The longitudinal velocity 



is 



vl = \\/3(a + QP) = \ y/3(C n + 2C 12 + 4C 66 ) . 



(127) 



Due to Propositions 15 and 16, a wide class of materials with the same S"-tensor and an 
arbitrary A-tensor will have exactly the same directions and the same velocities of the lon- 
gitudinal waves. The elasticity tensor for such materials can be written as 



a (3 + 2 Pl (3 + 2p 2 2p 3 

* a (3 + 2p 4 

* * a 

* * * (3 — P4 

* * * * 

* * * * 





2p 5 

2p 6 

P ~ P2 ~P3 

* P ~ Pi 



(128) 
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A Irreducible decomposition of tensors of rank p 

In order to understand what should be considered as the proper decomposition of a tensor, 
which is of rank p in its covariant or contravariant indices, that is, T l ^ kl or T^, we must look 
closer on its precise algebraic meaning. Let us start with a vector space V over the number 
field F. This construction comes together with the group of general linear transformations 
GL(dimV, F). This group includes all invertible square matrices of the size dimV x dimV 
with the entries in F. In our case, F — R and V = M 3 . Consequently, we are dealing with 
the group GL(3,M) that can be considered as a group of transformations between two bases 
{ei} and {e v } of V: 

e v = L/ ei V E GL(3, E) . (129) 

The vector space V cannot be decomposed invariantly into a sum of subspaces. Indeed, 
every two nonzero vectors of V can be transformed each other by the use of certain matrix 
of GL(3,M). Thus, the group GL(3,M) acts transitively on V. 

A tensor of rank p is defined as a multi-linear map from the Cartesian product of p copies 
of V into the field F, 

T : V x ■ ■ ■ x V -» F . (130) 

p 

The set of all tensors T of the rank p compose a vector space by itself, say T. The dimension 
of this tensor space is equal to n p . Thus, for rank 4 in 3d we have 3 4 = 81. As a basis in T, 
we can take tensor products of basis elements in V, 

e h ®e i2 ®---®e iv . (131) 

Accordingly, an arbitrary contravariant tensor of rank p can be expressed as 

T = T ili2 -^ e h ® e i2 (8) • • ■ ® e lp . (132) 

Under a transformation ( I129p of a basis of the vector space, the basis of the tensor space 
is multiplied by a product of the matrices This is a 'derived transformation' in the 

sense of Littlewood (1944) [21]. An important fact is that the tensor space T is decomposed 
to a direct sum of subspaces which are invariant under general linear transformations. For 
instance, the span of the basis 

e (il (g) e i2) (g) • • ■ <gi e ip (133) 

describes the subspace of all tensors which are symmetric under the permutation of the two 
first indices. Another subspace is obtained as the span of all basis tensors of the form 

e [h <g> e i2] (8) • • • <8) e ip , (134) 

which is the subspace of all tensors antisymmetric under the permutation of two first indices. 

By taking different permutations of the basis vectors we can obtain different subtensors 
(elements of the subspace of T) of a tensor of arbitrary rank. There are all together p\ 
permutations for the system of p objects. In fact, every GL(3,M) -invariant subtensor can 
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be obtained by the operator which permutes the indices, see Littlewood (1944) [21]. Since 
we need the sums of all the tensors with permuted indices, we must extend the group of 
permutations to its group algebra (Frobenius algebra). This construction involves formal 
sums, in an addition to the group multiplication. 

In this way, each tensor of rank two or greater can be irreducibly decomposed under the 
action of the linear group GL(3, R). The original tensor is written as a linear combination 
of simpler tensors of rank p, which, under the action of GL(3,M), transform only under 
themselves. These partial tensors obey their specific symmetries in addition to the symmetries 
of the initial tensor. Possible types of the irreducible tensors are determined by the use of 
Young's tableaux. 

A.l Example: Young's decomposition of second rank tensors 

Consider a generic second rank tensor which is decomposed into a sum of its symmetric 
and antisymmetric parts, 

rpij _ rp(ij) _|_ rp[ij] (135) 

It means that one starts with the tensor space 

T = Span {ej <S> ej} (136) 

The span of the linear combination of basis tensors 

S = Span{(e(j © e^)} (137) 

composes the subspace of tensors symmetric under the permutation of two indices. Indeed, 
an arbitrary tensor in S is decomposed as 

S = S ij e {l © . (138) 
Thus we have = 0. Another subspace A is obtained as a span of the linear combination 

A = Span{e[j <8> e^} . (139) 
This is a subspace of antisymmetric tensors, 

A = A ij e {i © e j} . (140) 

Here A^ = 0. In this way, the tensor space is represented as a direct sum of its subspaces 
T = S © A. These subspaces are invariant under the transformations of GL(3,R) while 
further decomposition into smaller subspaces is impossible. Hence we have an irreducible 
decomposition. 

In Young's description, this decomposition is given by two diagrams that are graphical 
representations of the permutation group S2 

n®n=m®H- m) 
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The left-hand side here denotes the tensor product of two vectors, i.e., a generic asymmetric 
second order tensor. The right-hand side is given as a sum of two second order tensors. The 
symmetric tensor is represented by the row diagram while the antisymmetric tensor is given 
by the column diagram. In our simplest example, only two tables given in (11411) are allowed. 
The next step is to fill in the tables with the different indices of the tensor. In particular, we 
write 

+ (ij))A ij = A ij + A ji , (142) 



1 



and 



Aij — A 



(143) 



Here / denotes the identity operator, (ij) is the permutation operator, (/+(r?)) and (/ — (ij)) 
are respectively Young's symmetrizer and antisymmetrizer operators. Thus, we have 



j 


= a 


i 


.1 



(144) 



The decomposition fl 1 3 5 [) is then obtained by inserting suitable leading coefficients a and (3. 
In general, these coefficients are calculated by using a combinatorial formula. It is easy to 
see that, for the completely symmetric and completely antisymmetric diagrams of n cells, 
the coefficients are equal to 1/n!. Thus, the coefficients in (1 144ft are equal to 1/2 and the 
decomposition H135[) is recovered. 



A. 2 Young's decomposition of fourth rank tensors 



The first step is to construct Young's tableaux corresponding to a fourth rank tensor. 

Rule 1. The four cells representing the indices of the tensor must be glued into tables of all 
possible shapes. The only restriction is that the number of cells in any row must be less or 
equal to the number of cells in the previous row. The number of irreducible subtensors of a 
given tensor is equal to the number of the tables of all possible shapes. 

Due to this rule, a generic fourth order tensor can be irreducibly decomposed into the sum 
of five independent parts. These parts are described by the following Young's diagrams, 
which are the graphical representation of the permutation group S4, see Boerner (1970) [1] or 
Hamermesh (1989) [Hj, 



(145) 



The left-hand side describes a generic fourth order tensor. On the right-hand side, the first 
diagram represents the completely symmetric tensor. The middle diagrams are for the tensors 
which are partially symmetric and partially antisymmetric. The last diagram represents a 
completely antisymmetric tensor. 

The next step of Young's procedure is to fill in the tables with the indices. In order to 
avoid the repetitions, the following rule is used: 
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Rule 2. In each row and each column of Young's table, the positions of the indices, i.e., 



the numbers k, i 2 



<*i>- 



are inserted in the increasing order. The symmetrization operators 



correspond to the rows of the table. Since different rows contain no common indices, the 
corresponding symmetrization operators commute. Different antisymmetrization operators 
correspond to the columns of the table. They also commute with each other and can be 
taken in an arbitrary order. The symmetrization and antisymmetrization operators act on 
the same indices, hence they do not commute and thus must be taken in some fixed order. 



Thus, due to graphical representation ( 1145P of the permutation group S4, we have the 
following parts of a generic three dimensional tensor T ili2i3i4 
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1 3 
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+ 



7 



k 


k 


k 


k 



+ 
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^l 



'2 



■'3 



«4 



(146) 



Rule 3. The coefficients a, /3, 7, 5, e are determined by the combinatorial formula, see Hamer- 
mesh (1989) [14J. They can be calculated using the Mathematica package 'Combinatorics 7 . 
The first and the last coefficients are especially simple. For an n-order tensor, 



a 



1 

n\ 



(147) 



For a generic 4th rank tensor, all these diagrams are relevant. In an n-dimensional space, 
it has n 4 independent components, which are distributed between the diagrams (I146p . The 
explicit form of the corresponding terms can be found in Wade (1941) [4"0]. 
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